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STABLE LONGITUDES FOR 12 HOUR 


ECCENTRIC ORBIT SATELLITES 

Carl A. Wagner 
Geod 30 iamics Branch 

Trajectory Analysis and Geodynamics Division 

ABSTRACT 

The accelerated longitude drift regimes of eccentric 12 hour orbits are con- 
sidered, due to the resonant geopotential. These quasi-stationaxy orbits are 
under investigation at the N.A.S.A. for large payload earth surveillance space- 
craft in the small applications technology satellite (SATS) program. Inclinations 
near 'critical' are especially attractive because of their possible long term sta- 
bility, and these are examined in detail. 'Stable' equator crossing longitudes 
for these satellites are found as functions of the argument of perigee. As long 
as the argument of perigee is between ±100^, these 'stable' longitudes lie within 
a narrow range. Maximum 'east-west' station keeping requirements for these 
satellites at non equilibrium positions are in the range of .3 to 3 m/sec/yr. 


Presented at: The First Goddard Astrodynamics and Geodynamics Conference; 

Greenbelt, Maryland, October 1971 
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INTRODUCTION 


Among the most attractive orbits for space applications is the eccentric 12 
hour orbit. This is a quasi-synchronous orbit which, with a low perigee, spends 
most of it's time within two small lattitude-longitude boxes at true synchronous 
altitudes (at about 36,000 km; see Figure 1 for example). Thus, if a mission 
requirement of an absolutely continuous link with a given ground station can be 
slightly relaxed, this orbit may be preferred over the more conventional 24 
hour 'geostationary' one. In fact the U.S.S.R. has already utilized this orbit ex- 
tensively for communications and other applications in their Molnisra and Cosmos 
series of satellites. Recently, NASA has considered utilizing such orbits in its 
small application technology satellite program. The first such 12 hour orbit, 
of Cosmos 41, was achieved in 1964. 

The payloads in 12 hour eccentric orbits are significantly heavier than in 
24 hour geostationary ones because (among other factors): 

1. The orbit is closer to the earth 

2. The orbit is highly eccentric 0.70) compared to the circular geostation- 
ary orbit, and 

3. No inclination change maneuver has to be performed for launches from 
high latitude sites (such as those in the Soviet Union). 

On the other hand, the circular 'geostationary' orbit is inherently more 
stable than the 12 hour-eccentric one whose perigee is subjected to strong long 
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term luni-solar perturbations. In fact maz^ of the Molniya's have already 

been lost (from rapid perigee collapse) due to improper initial placement of the 

[2] 

orbit node with respect to the moon's node. However, with proper node 
placement, this early demise can be avoided. Cosmos 41 is still in orbit in 
1971 and, according to Quinn et. al.,^^^ will stay in orbit for another 39 years. 

Radiation pressure also affects the eccentric orbit to a greater extent than 
the circular. However for the usual satellite payloads, where the projected area 
to mass ratio is considerably less than 1 cm^/gm, these perturbations do not 
affect the orbit substantially. 

The dominant perturbation affecting the longitude of the 12 hour orbit is due 
to resonance with longitude terms in the earth's gravity potential since the geo- 
graphic trace (latitude, longitude and altitude) ' 3 repetitive every day. Because 
of it's eccentricity, the resonant regime is more complex than for the 24 hour 
circular orbit. However it's stability is only affected by the significant rotation 
of the line of apsides, due to the earth's oblateness As a consequence of 
this rotation only two general 12-hour eccentric orbits exist which have long 
term longitude stability. These are the equatorial orbits, discussed in detail 
by Allan in 1967, and the critically inclined orbit (at 63.4** and 116.6°) 
briefly described by Wagner in 1968. lii both of these cases a stationary 
disturbing potential on the orbit can be constructed. This insures 
that, while the resonant libration regime is more complex than in the circular 
orbit case, the motion under this potential is always periodic. In particular, for 
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the equatorial and critically inclined orbit, cases of stable equilibrium can be 
found so that satellites placed in these orbits will not need stationkeeping to 
maintain a given geogp:apIdc trace* For the general incliimtion case, rotation cx 
the line of apsides gives rise to a non stationary potential, aperiodic motion 
pnd no stable longitude placement* This case has been discussed by Gedeon, 
Douglas and PEdmiter^^^ in 1967* It is well illustrated by the mean longitude 
history of the satellite 1966 96A (Intelsat 2- FI) shown in Figure 2* (The mean 
longitude (\) in this figure is the combination (M + cd )/2 + N - where M, 
and N are the satellites mean anomaly, argument of perigee and right ascen- 
sion of the ascending node, and is the earth's Greenwich hour angle*) For a 
circular 12 hour orbit this longitude would also be the actual geographic ascend- 
ing equator crossing longitude* For an eccentric 12 hour orbit there can be a 
considerable difference between this mean longitude and the actual equator 
crossing* (This difference will be described later as a function of the argument 
of perigee and the eccentricity*) Clearly the longitude excursions of this satel- 
lite are not periodic but are increasing with time as predicted by Gedeon et al.^^^ 
but contrary to the expectations of Allan*^^^ (Allan felt this orbit of 18 degrees 
inclination was sufficiently clo<^e to equatorial to have a stationary geopotential 
disturbing potential*) By contrast, the longitude history of 1964 49 F (Cosmos 41 
rocket) wifii an inclination near 'critical,' is nearly periodic (see Figure 3)* 

It's excursions are characteristic of a widely swinging pendulum* 
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Since 1969, starting with Molniya 11, the U.S.S.R. has placed their closely 
commensurate 12 hour satellites in orbits with ascending equat:>r crossings close 
to 112 degrees west. This longitude placement has proved to be stable for a 
wide range of perigee arguments. The shallowest excursions about this appar^t 
equilibrium position was shown by Molniya 12 in 1969-1970 before it's perigee 
collapsed in the fall of 1970 (Figure 4). 

This paper will derive all the resonant equilibrium longitudes for these 
satellites as hmctions of inclination and argument of perigee. Only the low 
perigee orbit (e = .725, perigee height = 930 km) is examined in detail because 
it is the most useful for the applications satellites. 
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ANALYSIS - THE DRIFT REGIME OF 
AN ECCENTRIC 12 HOUR SATELUTE 


The equation of motion of the mean longitude (\) of a closely commensurate 

r 1 f 1 

resonant satellite due to the geopotential is given by Wagner'^ ' as: 


k = 127T* 




rtlevant 









A 

0 


-cos m 1 A. + 1 

- 1 / - m even 


\ m / 


ma-^(FG)^„ 
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V m J 



radians/sidereal day^. 


( 1 ) 


In equation (1); the and are the unnormalized harmonic coefficients of 
degree i and order m > m) in the geopotential expansion/^^ a is the orbits 
semimajor axis in earth radii, and the (FG)^^ and 7 ^ are amplitudes and phase 
angles of composite harmonic vectors which depend on the inclination eccentricity 
and argument of perigee of the orbit. 


The composite vectors (FG) o are determined from Kaula*s F and G func< 
tions and his decomposition of the disturbing geopotential potential m) into 
{i , m, p, q) components The amplitudes (FG)^^ are given as; 
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FQp cos(-qoj) 


E 




p , q resonant 


p.q resonant 
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The phase angles are given as: 


n™ = T®" 


.1 


22 


p,q, reionant 


(?') 


i: 

p.q rtsonant 




The relevant resoiiant terms {i, m, p, q) are determined from the Indlcial 
equation: 


^ - 2p + q = m/s, ( 4 ) 

where s is the rational fraction representing the commensurability of the orbit 
in revolutions per day. For the 12 hour orbit, s = 2 and the dominant components 
on the highly eccentric satellite (e .7) are: (2, 2$ 1, 1), {2, 2, 0, -1), 

(3, 2, 1, 0), (4, 2, 1, -1) and (4, 2, 2, 1). Together these terms account for all 
but at most 10% of the resonant acceleration on the 12 hour satellite. All but 
at most 20% of the acceleration on this orbit is accounted for by just the two (2,2) 
components. Furthermore, for orbits in the vicinity of critical inclination all 
but about 25% of the acceleration is accounted for by the single term (2, 2, 1, 1). 
Similarly, for equatorial orbits, the single term (2, 2, 0, -1) accounts tbr all but 
about 25% of the acceleration. Therefore it is both instructive and useful to 
work out these two simplest cases first, for low and high inclinations, as fairly 
good approximations to the drift regpime. 
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The analysis that follows will only be concerned with determining the con- 
ditions for longitude equilibrium on these satellites* It must be bom in mind 
that only fo;' the cases near 'critical' inclination (33.4*) and for equatorial orbits, 
can £ach equilibrium conditions be maintained for long periods of time. As il- 
lustrated in figures 2 and 3, an 18* inclined orbit does not maintain it's stability 
in this sense while a 68 degree orbit does. 

Low Inclination Approximation 

In this case, the equation of motion, (1), of the longitude can be reduced to the 
simple pendulum equation (for constant ^ , a, 1 and e), involvixig the term (2, 2, 0, 
-1) only: 


K = [Cj sin (2k + ci;) + Cj cos (2k + co)] . ( 5 ) 

The zero's of this equation give the equilibrium conditions for the orbit. Equili- 
brium exists when: 


j - [C j sin cc) + Cj cos o) ] 

\ Tan" — T ; 

2 C j cos CO - Cj sin CO 


But Cj = KC 2,2 ^2 = ~KS 2 2 f where; K = 2477^ ^ (I) ^ (e)/a, for the 

12 hour orbit. Also, we have from the definition of the harmonic coefficients 

^ 2 , 2 "" *^ 2,2 ^^ 2 . 2 * ® 2,2 *^ 2,2 ^^ 2 , 2 » whero, from current values, 

^ 2 , 2 ^ With these substitutions for and C 2 , the equilibrium longitudes 
are given as: 
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( 6 ) 


1 ^ 1''°® sin a. 

2 L^os 2^22 

1 ^ .t si" (2^22 -") 

= ~ Tan ' 

2 cos (2^22 " 

Equation (6) has four solutions: 

= >^22 -^/2, ^22 + 90®, ^22 + 180® and ^22 “^/2-< 270®. 

However since the 12 hour orbit always has two mean longitudes 180° apart 
(corresponding to the two ascending equator crossings), only two of these solutions 
are distinct, for example: 


- sin 2 \j 2 
+ sin 2A-22 ® 


os co” 
i n 


= ^2, 2 ” and ^2, 2 ” + 90 ® . 


( 7 ) 


One of these equilibrium orbits is stable and the other unstable. To find which 

90 00 

is which, we note (from Figure 5) that the zeros of X where d Vd^- < 0 are stalTe 
longitudes, while those for which 6K/dk > 0 are unstable. Differentiating (5) 
with respect to k and substituting the first i^olution of (7), we have: 


dk 


= 2[Cj 


COs 2 \- 2,2 


-C2sin2\2,2] =2KJ2.2 


[ CO S^ 2 X 2,2 + 2^2,2^ 


= 2KJ 


2 , 2 * 



since J > 0, d\/dX for this first equilibrium solution has the sign of K, or 
• • 2 

the sign of Fj 2 . 0 ^ 2 . 2.0 (I) = 3 (1 + cos I)^ /4 :> 0, and 

G 2 . 0.-1 (©) = + e^/16 + . . . ^ 0 (see Reference 9). Therefore, “ ^ 2.2 

is a stable equilibrium longitude and similarly it can be shown that 
\ . -oj/2 + 90^ is an unstable equilibrium longitude. 

R. R. Allan in 1967 placed the stable perigee longitudes for the equatorial 
12 hour eccentric orbit at ^2 2 ^ 22 *^ redundant position). 

We can check this placement by the calculation above with reference to 
Figure 6. For an equatorial orbit, the geographic longitude of perigee is 

I 

77^ + N- ^g, when the satellite (s) is at perigee (or M = 0). But when the 

satellite is at perigee, the mean longitude is = Ci^/2 + N - This longitude 
is stable when \ = ^2 2 ' Therefore, the stable condition implies: 

^2,2 “ ^ 2 -I- N - ( 9 ^ ; or 

^2 . = " + N-e^ = 7T^ , 

as first postulated by AUan.^^^ 

High Inclinat on Approximation 

In this case, the equation of motion of the mean longitude, (1), involves only 
the term (2, 2, 1, 1) and is formally the same as (5) with >a> substituted for a; in 
the sine and cosine arguments. Equilibrium values are therefore at: 
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^ = ^2.2 + and X. 2,2 + ^/2 + 90 ° . 


( 8 ) 


By a calculus entirely similar to the previous case, recognizing that: F 2 2 , i (1) 
= 3 sin^ 1/2 0, and 62 ^1 (e) = 3e/2 + 27 e^/16 + 0, it can be shown that 

= ^ 2.2 + ‘^/2is the unstable equilibrium longitude and e “ ^ 2 . 2 + ^ /2 + 90° 
is the stable equilibrium longitude. Figure 7 shows these approximate locations 
of the stable mean longitudes for low and high inclinations. 

For intermediate inclinations the situation is more complicated. As 
mentioned before no geostationary disturbing potential exists for the general 
inclination case. Therefore the equilibrium longitudes derivable for them are 
not valid over the long term. In fact the time of validity depends on ^ amount 
of rotation of the line of apsides. This can be appreciated from Figure 7 since 
for each of these harmonic terms, a full range of perigee arguments implies a 
full range of 'stable' longitudes. For intermediate inclination 12 hour orbits, 
the period of apsidal rotation is of the order of a few years (Figure 8). In ad- 
dition,for intermediate inclinations, more than one harmonic term has appreci- 
able influence. But for inclinations fairly near critical, which are most useful 
in the applications, only the two terms (2, 2, 0, -1) and (2, 2, 1, 1) carry all but 
about 10% of the effect. Therefore it will be instructive and practical to calcu- 
late the approximate quasi-equilibrium longitudes from this particular 
combination. 
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Intermediate Inclination Approximation 


The equation of motion (1) including the two terms (2» 2, 0, -1) and (2, 2, 1, 1) 
is: 


\ ^ (FG)j_j [Cj j sin(2A. + 7 j,,) -S, jCOs(2\ + 7j _ j)] , 


radians/sidereal day^, 


( 9 ) 


where (FG >2 j ^ 2.2 functions of I, e, and The vector interpretation 
of these functions [defined in equations (2) and (3)] permits them to be easily 
resolved by trigonometry (Figure 9). The individual component amplitudes 
(FG's) in this figure all are absolute values. Thus, using the law of cosines we 


have: 


(FG),,, = + [FG 2 _ 3 _,_j+FG^_j, .,- 2FG,_,_j.,FG, cos..] 


1/2 


( 10 ) 


Using the law of sines: 


FG< 


2 . L . ? . , determines 7, 2 

sin (-722"“^^ i:^n 2 aj 


^2,2 = "i 


f • P 

CO + sin 


*'G2.2,o..i Sin 20) 




(FG) 


2,2 


( 11 ) 


I .-If sin 2 CO 

= - Jco + sin 1 

[ [(X2 + 1 - 2Xc< 


os 2cjS) 


whereX= IFG 2 2 ,i.iI/IfG 2 , 2 ,o,-J • Table 1 gives values of x as a function of 
inclination taken from the graphs in reference 10. 
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The zero's of (9), giving the 'equilibrium' longitudes, are found from the 
equation: 

cos 2^22 sin( 2 \ + ^2,2) = sin 2^2,2^°®^^^ ^ ^2,2^’ 
or; 

Tan (2 \ + 72 , 2 ^ ~ ^^2, 2 ’ 

The 2 unique solutions to this equation are: 

^ 2, 2 " ^2, 2^^ ’ ^ 2, 2 " ^2, 2^^ ^ 

It can be shown that the second solution is the one for stable equilibrium. 
Thus, for the intermediate inclination case (second approximation): 

(Stable) =^2,2 - S(o;,I), ( 13 ) 

where S is one half the arc sine function in equation (11), Ms a small sinusoidal 
function for X >> 1. For example, for the most useful case of the critical inclin- 
ation satellite of high eccentricity (e = .725, perigee height = 930 km); X = 6,20 
and § (^, 63.4°) is plotted in Figure 10. 

Figure 11 shows the variation of (stable) with argument of perigee, for 
intermediate inclinations, under various approximations. The simplest, straight 
line solution is most accurate for high inclinations and considers the (2, 2, 1, 1) 
term only. The second approximation (including the (2, 2, 0, -1) term) is shown 
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for the critical inclination case only. Third approximations (including the next 3 
most influential terms (3, 2 , 1, 0), (4, 2 , 1, -1) and (4, 2, 2, 1)) for I = 50% 60“ 
and 70® are also shown in Figure 11. These are calculated directly from the 
zero’s of (1) for these terms. The gravity field used for these calculations is 
the Smithsonian Standard Earth I, Ml Field.^''^ Even at 50® inclination the 
maximum deviation from the straight line solution is only 17®. Consideration of 
higher degree ano order terms will add or subtract less than 10% of 17® or less 
than 2® to those equilibrium longitudes. Inclinations less than 50® are not con- 
sidered in Figure 11 because, as mentioned previously, the apsidal rotation is 
too rapid for the ’equilibrium’ longitudes to be considered anything but momen- 
tary values. However, the near* equatorial satellite is a special case, as Allan 
found. For it, the apsidal rotation is in the same sense as the earth’s rotation 
and a stationary geographic trace is assured with proper adjustment of the mean 
motion. As inferred previously, near equatorial 12 hour eccentric orbits are 
stable if one of their periapses is within a few degrees of x ^ ^ or 15® west. 

For the applications it remains to calculate the actually geographic position 
of these ’stable’ orbits when the ’stable’ mean longitudes are given (as in Figure 
11). This will be done in the following section by finding the relation between 
the actual ascending equator crossing (AEC) longitude and the mean longitude (^ ) 
for a commensurate (or resonant) orbit. 
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THE ASCENDING EQUATOR CROSSING LONGITUDE 
Consider the arguments in the orbit and equatorial planes from Figure 12. 
y is the vernal equinox, G is the position of Greenwich at epoch, when the satel- 
lite is at SAT, Go is the Greenwich position at a time Tp + T^ earlier, when 
the satellite is at the ascending node (AN), is the right ascension of Green- 
wich, ^ is the argument of perigee and M is mean anomaly, all at epoch. When 
the satellite is at AN, is the longitude of the ascending equator crossing. 
We make two simplifying assumptions; one, that over the time Tp + T^ the 
node and perigee are stationary and two, the commensurate mean motion of the 
satellite, n, is s^.. 

From Figure 12: 


(Tp +Tc). 

But since Tp = M/n, and n = s 

^AEC=N-^e + l/s(M + T^n). (13) 

The quantity T^n may be called the mean anomaly equivalent r,l o;,(a;j^), (i.e. it 
is the angle swept out at the average rotation rate in the actual time for the 
satellite to go from AN to tt.) We recall^®^ ’ that the definition of the mean 
longitude (A.) for a general commensurate orbit is: 

\ = + 1/s (M + o;) (14) 
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Comparing (13) and (14): 


^AEc = ^ + 1 /s (^M - (15) 

Let us call the quantity ^ ^ perigee exces^^ (Aoj). Evidently Ao; is zero 
for a circular orbit and also zero for eccentric orbits when ^ = 0^ and 180*. 

A a; is clearly some sinusoidal function of cj for a given e. The function Ao) can 
be calculated from Keplers equation: 

M = E - e, sin E , 

and the relationship between the true anomaly (f) and the eccentric anomaly (E), 

E= 2 Tan-‘ [(5^)‘"Tan f/ 2 ] . 

In terms of the mean and eccentric anomaly equivalents of perigee (oj^ anda;^), 
these relations become: 

co^ = a)g - e sin and (17) 

Aoj =: 0^1^ - CO . (18) 


Tan co /2 


I 

- 


(16) 
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It is noted that A a; is an *odd' function about a; = 180**. The asymptote of A a;, 
when 6 = 1, are the parallel lines of -1 slope through a; = 0* and 360*". For use 
in applications fore<<l,Ao; = -e sin . For any eccentric commensurate 
orbit, the graph of Aoj as a function of oj and e is given (from evaluation of (16), 
(17) and (18)) in Figure 13. 

Finally, figures 11 and 13 are combined (with equation (15)) to give (in 
figure 14) the 'stable* equilibrium equator crossing longitudes on 12 hour ec- 
centric orbits. These are given for a range of inclinations useful for applications 
satellites (with an eccentricity = .725, giving a perigee altitude essentially out 
of the atmosphere at 930 km.) In Figure 14 is also shown the observed and 
computed 'stable* equator crossing longitudes for Molniya 12 which confirms the 
validity of this chart to within a few degrees. 
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STATION KEEPING AT NON EQUILIBRIUM LONGITUDES 


Figure 15 shows the longitude acceleration on the maximum altitude 12 
hour orbit of 50* inclination for = 270*, a typical orbit that has been considered 
advantageous in the SATS program/ Only the dominant terms '(2, 2, 0, -1), 

(2, 2, 1, 1), (3, 2, 1, 0), (4, 2, 1, -*1) and (4, 2, 2, 1)) have been included in the 
evaluation of the acceleration from equation (1). The Smithsonian Standard 
Earth I, Ml gravity field has been used in this computation. 

It is interesting that the maximum resonant acceleration on this 12 hour 
satellite is of the same order of magnitude (1 to 10 x IC"^ rad/day ) as that on 
the 24 hour satellite.*- This rough *equipartition* v f the disturbing force among 
all the *one day* resonances was first noted by Douglas et al. in 1968.*-^®^ *One 
day resonances have ground traces which repeat in one day. In reference [12] , 
the maximum station keeping requirements on 24 hour satellites were estimated 
to be about 2 m/sec/yr. Because of equipartition we can expect the requirements 
to be the same order of magn^^uie for all one day resonances. However, some 
fuel savings can be realized . =)ccentric orbits by performing the maneuvers 
at perigee where a given change in velocity (AV) has the greatest effect on the 
energy of the orbit. The following estimate of the requirements assumes that 
all station-keeping thrusting occurs at perigee. 

If we ask *what is the longitude acceleration (x) in a commensurate orbit 

due to a continual change in orbit period (T) of A T each commensurate period 

[u] 

of n days the answer is exactly the same as for the circular case: 


IV 



AT 2 

\ = - 2tt n * — radi ans/(n ' days) 


= - 2t 7 (1/n ') 


AT 

— radians/day^ 


This equation could be written without specified dimensions with & ^ - 2 ' tt , But 
from Keplers 3rd law; 


AT _ 3 Aa 
T "" 2 a ’ 

so that 

K = ^ Aa rad/day^. 

a 

In the above equation the negative sign means that an increase in the semimajor 
axis causes the orbit to fall behind the rotating earth. We now need to determine 
only the AV required to change the seminmjor axis A a to complete the station 
keeping specifications. For uds purpose we use the two-body energy integral 
for elliptic motion in the vis-viva form: 

where ^ is the central gravitational constant (3.986 x 10^ km^/sec^), r is the 
distance to the sateUite from the center of the earth and V is the satellite's 
velocity. Taking differentials of (19) with respect to V and a (holding r constant, 
at perigee (p))$ 
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2V AV 

p p 


/i. A a 


a 


2 


But from (19) again (with Tp = a (1 - e)); 




Therefore, perigee velocity change (^Vp) results in a semimajor axis change, 
A a, according to the formula; 


Aa = 2(a^//x) 


1/2 ft + e\ 


1/2 


AV . 

p 


Substituting this equation into the acceleration equation (and reversing the sign) 
gives the station keeping requirements (for perigee maneuvers) as: 


AV = 

p 


(n') (M/a) 


6tt 




1/2 


K . 


(20) 


Equation (20) gives the rxiinimum A V requirements for a commensurate period 

t* 

of il' sidereal days in length/time units, if X is given in units of radians/ 

2 , 

sidereal day and n is given in sidereal days. The sign is reversed because 

the station keeping maneuver must be in an opposite sense to the A V rate 

•• 

equivalent to the acceleration, X , caused by the geopotential. Therefore, the 
yearly minimum Ay requirement for station keeping is 366 /n' times A Vp or: 
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AV (m i n i mum ) = 


366 

6tt 



(21) 


= 1.225X 10^ 


1 - 

.1 + Si 


1/2 •• 

m/sec/y r. , 


If a is given in km and K in radians/sidereal day^. 

For example, the maximum altitude 12 hour orbit has the specifications 
a =■ 26,550 km, and e = .725. If X. = 5 x 10‘® rad/days, AV (minimum), (from 
(21)) is found to be 1.51 m/sec/yr or 4.95 ft/sec/yr. This calculation establishes 
the scale of the velocity requirement in Figure 15. Positive station keeping Av*s 
are to be made along the track, while negative Ay^s are to be made against the 
track of the orbit. 


As another example, consider the orbit proposed in Figure 1 for the small 
applications technology satellite (SATS) program. This orbit (of 50® inclination) 
has an argument of perigee of about 270® and ascending equator crossing (AEG) 
longitudes at 136® west and 45® east. 

From Figure 14, the 'stable* ascending equator crossing (AEG) longitudes 
for this (maximum altitude) orbit are at 60® east and 120® west. Thus the mean 
longitudes (^) of this orbit are displaced 15® west of the 'stable' mean longitudes. 
From Figure 15 (or Figure 11) the 'stable' mean longitudes of this orbit are at 
27® and 207® east. Therefore, the mean longitudes of the given orbit are at 
27 - 15 = 12® east and 207 - 15 = 192® east. 
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A simpler calculation of the mean lon^Tltude ( k) of this orbit can be made 


with the use of equation (15): 


^ * ^AEC 


1 

S 


Ao;. 


For this 12 hour orbit; s = 2, e = 0,725, ^ = 270® and (from Figure 13) = 66®, 
Thus, \ = 45 - 66/2 = 12® or 192® as found previously. From Figure 15, the 

—5 2 

resonant acceleration on this orbit is 5 x 10* rad./day and the station keeping 
requirements are 1.5 m/sec/yr. 
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SUMMARY AND CONCLUSIONS 


The eccentric 12 hour orbit has already been one of the most useful orbits 
for earth oriented applications. The United States Intends to make further use 
of these orbits for many applications in the decade of the 1970's and beyond. 
Since these orbits are strongly perturbed by resonant terms in the geopotential, 
it is of interest to estimate these effects, and, if possible minimize them for a 
particular mission. 

In general it is found that: 

1. The longitude acceleration on these quasi- stationary satellites is highly 
dependent on the argument of perigee. 

2. For each (equatorial and critically inclined) orbit specification (e, a;) 
only two longitudes for the placements of the equator crossing result in 
orbits which are in equilibrium with the resonant geopotential forces. 

One of these orbits is in stable equilibrium with these forces and the 
other is in unstable equilibrium with them. 

3. For intermediate inclinations, the 'sta de* longitudes change slowly with 
time (as the argument of perigee changes) as long as the orbit inclination 
is within about 15 degrees of critical. The stability of the equatorial 
orbit rapidly degenerates for inclinations above 10 degrees. 

Specifically it is found that for the maximum altitude 12 hour orbit of inter- 
mediate inclination (50® -» 70®): 
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1. 'Stable* orbits exist for ascending equator crossing longitudes between 
55** and 85** east, as long as the argument of perigee is within dilOO*. 

2. The majdmum station keeping requirements to keep the geographic trace 
constant against resonant geopotential forces are in the range of .3 to 3 
m/sec/yr. 
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Table 1 


Harmonic Amplitude Ratio (X) as a Function of Inclination 
for High Altitude (High e) 12 Hour Satellites 

I X Data Derived from 

0® 0 Referenco 10. 

10 0.123 

20 0.509 

30 1.16 

40 2.13 

50 3.54 

60 5.34 

(Crit. Ihcl.) 63.4® 6.20 

70® 9.90 
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1. Alternate Ground Trace for 12-Hour Orbit - 50® Inclination 

2. Mean Longitude (\) History for 1966 96 A (Intelsat 2 FI) 

3. Mean Longitude (X) History for 1964 49E (Cosnos 41 Rocket) 

4. Actual Ascending Crossing Longitudes (^aec) 

(Molniya 12) 

5. Stability Conditions With Respect to the Longitude (\) Drift 

6. Equatorial Orbit Reference Points 

7. Stable Equilibrium Longitudes (\) for High Altitude 12 Hour Satellites 

8. Period of Apsidal Rotation (T ) for a High Eccentricity 12 Hour 

p 

SateUite 

9. Composite Harmonic Functions for Intermediate Inclination 
12 Hour Satellites 

10. Equilibrium Mean Longitude Change as a Function of Argument of 
Perigee for a 12 Hour Satellite 

11. Stable Equilibrium Mean Longitudes (\_) for Maximum Altitude 12 
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Hour Satellites of Intermediate Inclination 


Figure 12. Reference Points of Orbit and Equatorial Planes on the Celestial 
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Perigee Excess (Aa>) for Eccentric Orbit Satellites 


^Stable' Ascending Equator Crossing Longitudes (^aec) Hour 
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Acceleration and Station- Keeping Requirements on a 12 Hour 
Satellite 
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Figure 2. Mean Longituoe (A) History for 1966 96 A (Intel sot 2 FI) 
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Figure 3. Mean Longitude (A) History for 1964 49E (Cosmos 41 Rocket) 
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Figure 4. Actual Ascending Crossing Longitudes for 1969 61A (Molniya 12) 
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Figure 5. Stability Conditions With Respect to the Longitude (A) Drift 
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Figure 6. Equatorial Orbit Reference Points 
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STABLE EQUILIBRIUM 
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Figure 7. Stable Equilibrium Longitudes (X) for High Altitude 12 Hour Satellites 
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Figure 8. Period of Apsidol Rotation (1^) for a High Eccentricity 12 Hour Satellite 
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Figure 9. Composite Harmonic Functions for Intermediate Inclination 12 Hour Satellites 
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ORBIT PARAMETERS 
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Figure 10. Equilibrium Mean Longitude Change os a Function of Argument of Perigee for a 1? Hour Satellite 
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Figure \ 2 . Perigee Excess (Acu) for Eccentric Orbit Satellites 
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